Refined similarity hypotheses in shell models of turbulence 
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A major challenge in turbulence research is to understand from first principles the origin of anoma- 
lous scaling of the velocity fluctuations in high-Reynolds-number turbulent flows. One important 
idea was proposed by Kolmogorov [J. Fluid Mech. 13, 82 (1962)], which attributes the anomaly 
to the variations of the locally averaged energy dissipation rate. Kraichnan later pointed out [J. 
Fluid Mech. 62, 305 (1973)] that the locally averaged energy dissipation rate is not an inertial- 
range quantity and a proper inertial-range quantity would be the local energy transfer rate. As a 
result, Kraichnan's idea attributes the anomaly to the variations of the local energy transfer rate. 
These ideas, generally known as refined similarity hypotheses, can also be extended to study the 
anomalous scaling of fluctuations of an active scalar, like the temperature in turbulent convection. 
In this paper, we examine the validity of these refined similarity hypotheses and their extensions to 
an active scalar in shell models of turbulence. We find that Kraichnan's refined similarity hypothesis 
and its extension are valid. 

PACS numbers: 47.27.-i, 47.27.eb 



I. INTRODUCTION 

Much effort in turbulence research has been devoted to 
the study of the possible universal statistics of velocity 
fluctuations in turbulent flows at high Reynolds number. 
The seminal vifork of Kolmogorov in 1941 (K41) J_] intro- 
duced the idea of universal homogeneous and isotropic 
statistics in high- Reynolds-number turbulent flows and in 
particular deduced that the statistical moments of the ve- 
locity differences, Sur{x,t) = \u{x + r,t) —u(x,t)\, where 
u is the velocity field, have power-law scaling with the 
separating distance r = |f]: 



{dur{x,ty-) ^ i{e)r) 



i/3 



(1) 



when r is within the inertial range. Here (. . .) denotes an 
ensemble average. The inertial range refers to the range 
of length scales that are smaller than those of energy 
input and larger than those affected directly by molecular 
dissipation, and the energy dissipation rate per unit mass 
e is given by 
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where v is the kinematic viscosity, x — {xi,X2,X3), and 
u — {ui,U2,u^)- The direct proportionality of the scal- 
ing exponents 7i/3 with the order n of the statistical mo- 
ments is equivalent to an independence of r of the func- 
tional form of the probability density function of 5ur. 
Experimental observations confirm the power-law scal- 
ing but indicate that the scaling exponents depend on 
the order n in a nonlinear fashion. This deviation of the 
velocity scaling behavior from the K41 results is known 
as anomalous scaling and implies that turbulent velocity 
fluctuations have scale-dependent statistics and are thus 
intermittent. 



A major challenge is to understand, from flrst prin- 
ciples, the origin of anomalous scaling. In his reflned 
theory in 1962 [2], Kolmogorov replaced the global av- 
erage energy disipation rate (e) by the locally averaged 
energy dissipation rate e^, given by 



r{x,t) 



Airr^ 



\h\<r 



e{x + h, t)dh 



(3) 



The possible dependence of (e" ) on r allows for a cor- 
rection to the n/3 scaling (see Sec. HII- As a result, this 
idea of Kolmogorov, which we refer to as Kolmogorov's 
reflned similarity hypothesis (RSH) attributes the origin 
of anomalous scaling to the variations of the local energy 
dissipation rate. Kraichnan Q later pointed out that the 
local energy dissipation rate e^ is not an inertial-range 
quantity and a proper inertial-range quantity would be 
the local energy transfer rate, Hr{x,t), which measures 
the rate of energy transfer into scales smaller than r 
at a local point x in space. In other words, Kraich- 
nan's idea, which we refer to as Kraichnan's RSH, at- 
tributes the origin of anomalous scaling to the variations 
of the local energy transfer rate. In statistically steady 
state, for r in the inertial range, (H^) = (e) = (cr) but 
Hr(af, t) ^ €r{x,t). Hence the two RSHs of Kolmogorov 
and Kraichnan would give different results for the anoma- 
lous scaling exponents of velocity structure functions. 

Anomalous scaling behavior has also been observed in 
the statistics of a scalar held advected by a turbulent ve- 
locity fleld. A passive scalar leaves the velocity statistics 
intact while an active scalar couples with the velocity 
and influences its statistics. The nonlinear problem of 
anomalous scaling of active scalars, like that of velocity, 
remains unsolved. A common example of an active scalar 
is temperature in turbulent convection in which temper- 
ature variations drive the flow. Both Kolmogorov and 
Kraichnan's RSHs can be extended to turbulent convec- 
tion to account for the anomalous scaling behavior of an 



active scalar. 

In this paper, we examine the vahdity of Kolmogorov 
and Kraichnan's RSHs and their extensions to turbulent 
convection in shell models of turbulence. Shell mod- 
els for homogeneous and isotropic turbulence have been 
proved to be very successful in reproducing many of 
the statistical features of turbulent flows, and in which 
high Reynolds number can be achieved with relative 
ease Q. Shell models for homogeneous turbulent con- 
vection have also been studied. We use two shell models: 
the Sabra model for homogeneous and isotropic turbu- 
lence [5|, which is an improved version of the so-called 
Gledzer-Ohkitani-Yamada (GOY) model [0,01, and the 
Brandenburg model for homogeneous turbulent convec- 
tion fs^. 

The rest of this paper is organized as follow. We first 
summarize the mathematical formulations of the RSHs of 
Kolmogorov and Kraichnan and the previous studies that 
examined their validity in Sec. [TTl In Sec. IIIIl we discuss 
how these RSHs are extended to turbulent convection in 
which temperature is an active scalar. We describe the 
Sabra shell model for homogeneous and isotropic turbu- 
lence and the Brandenburg shell model for homogeneous 
turbulent convection in Sec.lIVIand Sec. El respectively. 
In Sec. IVH we examine the validity of the RSHs of Kol- 
mogorov and Kraichnan for homogeneous and isotropic 
turbulence in the Sabra model and the validity of their 
extensions to turbulent convection in the Brandenburg 
model. Finally, we conclude in Sec. IVIII 



II. KOLMOGOROV AND KRAICHNAN'S 
REFINED SIMILARITY HYPOTHESES 



Kolmogorov's RSH can be stated as: 
5ur{x,t) = (j) [er{x,t)rY/^ 



(4) 



where is a random variable independent of r and statis- 
tically independent of e^ when Reynolds number is much 
larger than 1. Equation ^ implies that 



{{SurY 



gV3^^n/3 



(5) 



For homogeneous flows, {erix,t)) — {e{x,t)) and is inde- 
pendent of r and x. For n ^ 3, (e" ) generally depends 
on r and this r-depcndence allows for a correction to the 
K41 scaling. Thus Kolmogorov's RSH attributes the ori- 
gin of anomalous scaling to the variations of the local 
energy dissipation rate. 

Similarly, Kraichnan's RSH can be stated as 



SUrix,t) ^ijj [\nrix, t)\r]^/^ 



(6) 



where i/; is a random variable independent of r and sta- 
tistically independent of H^. The local energy transfer 
rate H^ can be defined using banded Fourier series [3|] 
Equation ([6]) then implies 



and Kraichnan's RSH thus attributes the origin of 
anomalous scaling to the variations of the local en- 
ergy transfer rate. As the scaling behavior of (|H,. |"/'^) 

would be generally different from that of (e" ), the two 
RSHs, Eqs.Q and ^, would give different results for 
the anomalous scaling exponents of the velocity struc- 
ture functions. 

Kolmgorov's RSH has been widely used in the dis- 
cussions of anomalous scaling of the velocity structure 
functions. There have been quite a number of exper- 
imental and numerical studies that examine the valid- 
ity of Eq. ([4]). Most of these studies checked whether 
6ur and e^ are correlated. Statistical correlation between 
6ur and e^ or its one-dimensional surrogate which repre- 
sents e by vidux/dx)"^ or [v/2){duy/dxY was indeed re- 
ported [alia, [ill, [il,[il- it was, however, noted ^^ 
that at least part of such a statistical correlation re- 
sults from kinematical constraints independent of Navier- 
Stokes dynamics. Moreover, it would be more direct to 
check the implication of Eq. ((3]) that the conditional ve- 
locity structure functions at fixed values of the e^ would 
be given by: 



{{5Ur 



^n/3gn/3 



(8) 



{{SurT) ^ (|n,r/3)r"/3 



(7) 



Calculations of such conditional velocity structure func- 
tions (using the longitudinal velocity difference instead 
of the whole velocity difference) were carried out using 
direct numerical simulations of isotropic turbulence jl6| 
but these simulations were limited to moderate Reynolds 
numbers making it difficult to draw definitive conclu- 
sions. In particular, a clear demonstration of a scaling 
behavior of r"/'^ at fixed e^ or a dependence of e" at 
fixed r is lacking. 

On the other hand, the validity of Kraichnan's RSH is 
much less examined. In a high resolution direct numerical 
simulation of isotropic turbulence [17| , it was shown that 
the scaling exponents of (IrHrl*'/'^) are close to those of 
the p-order transverse velocity structure functions. Since 
the transverse velocity difference is more intermittent 
than the longitudinal velocity difference [l^, this result 
implies that the scaling exponents of (|rHj.|P'''^) are close 
to those of the p-order whole velocity structure functions 
and is thus consistent with Eq. ([7]). However, a direct 
examination of Eq. ^ by studying the behavior of the 
conditional velocity structure functions at fixed values of 
the local energy transfer rate is yet to be performed. 



III. REFINED SIMILARITY HYPOTHESES 
FOR TURBULENT CONVECTION 

In turbulent convection, the statistics of 5ur as well as 
those of the temperature difference 5Tr = T{x + r,t) — 
T{x,t), where T{x,t) is the temperature field, are of in- 
terest. Buoyancy is expected to drive the dynamics, and 
thus affecting the statistics of 6ur and STr, for length 
scales r greater than the Bolgiano scale [13, llOl • When 



buoyancy is dominant, it was suggested [2l| that the 
statistics are governed by a cascade of temperature vari- 
ance, which isproportional to entropy in Boussinesq ap- 
proximation |22l | . Thus the role of the energy dissipation 
rate e is now played by the temperature (variance) dissi- 
pation rate or the entropy dissipation rate x, given by 



X 






(9) 



where k, is the thermal diffusivity of the fluid. In the same 
spirit of deriving Eq. ([¥]), 5uj. and 5T^ are expressed as 
functions of ag, Xr and r only, where a is the volume 
expansion coefficient of the fluid, g is the acceleration 
due to gravity, and Xr is the locally averaged entropy 
dissipation rate, given by 



Xr{x,t) 



A-Kr^ 



x{x + h, t)dh 



(10) 



h\<r 



As usual, the functional dependence is obtained by di- 
mensional analysis and the results are: 

5ur = KiagflXrixM^r^" (H) 

5Tr = ^T[agr^"'[xr{xM"'r^"' (12) 

Here $„ and $7 a-re random variables independent 
of r and statistically independent of Xr- Equa- 
tions (fTTjl and (I12p are the extensions [2^ of Kol- 
mogorov's RSH [Eq. (j4|)] to turbulent convection where 
buoyancy is dominant, and they attribute the origin of 
the anomalous scaling to the variations of the local en- 
tropy dissipation rate [24| . In particular, the conditional 
velocity and temperature structure functions at fixed val- 
ues of Xr would have simple scaling behavior in r: 



{{SUrT I Xr) 
{\STrr I Xr) 



r."/5 



(13) 
(14) 



given by that of Bolgiano-Obukhov (BO) [ig|,l20|. 

To get the corresponding expressions for the extensions 
of Kraichnan's RSH [Eq. ^] to turbulent convection, 
one replaces Xr by the locally entropy transfer rate Hf, 
which is defined as the rate of entropy transfer into scales 
smaller than r at a local point x in space. The results 
thus read: 

Sur = #„(a5)'/'|n^(a?,t)|'/'r3/5 (15) 

STr = ^T(a5)"'/^|nf(f,t)p/Vi/5 (16) 

Here XPu and ^t sue random variables independent of r 
and statistically independent of H^ . Equations ^T5\i and 
([TS]) therefore attribute the origin of anomalous scaling 
to variations of the local entropy transfer rate. Their 
implications are that the conditional velocity and tem- 
perature structure functions at fixed values of H^ would 
have simple BO scaling: 



{{SUr)'" I H^) 

{\STr\" I n^) 



^3n/5 



(17) 
(18) 



Turbulent convection is often investigated experimen- 
tally in Rayleigh-Benard convection cells heated from be- 
low and cooled on top (see, e.g., [2^, [2y, [23 for a re- 
view). Such confined convective flows are highly inho- 
mogeneous in which thermal and viscous boundary lay- 
ers are present near the top and the bottom of the cell, 
and coherent structures are present. It was argued [23| 
that the presence of buoyant flow structures known as 
plumes could affect the scaling behavior, causing BO 
scaling to be invalid. Indeed there were indications from 
direct numerical simulations J28| and analyses of experi- 
mental data [29] that the scaling behavior of the central 
region of confined turbulent convection is not well de- 
scribed by BO scaling plus corrections. We note that in 
Ref . [2^ , the statistics of Sut = \ u{x, t + r) — u{x,t)\ and 
STr = T{x,t -\- t) — T{x,t), the temporal counterparts 
oi Sur and 5Tr, were studied because only measurements 
taken as a function of time at a fixed point in space were 
available. The validity of Eq. (|12p was explored [2J| by 
studying the conditional statistics of 6Tr at fixed val- 
ues of Xr, estimated by Xt cx (1/t)/^ '^ n{dT/dt')'^dt' . 
It was found that for scales larger than the Bolgiano 
scale, {\5TrY'\xr)l{{\Tr?\XT))'^''^ become independent of 
T. This result indicates that the scale- or r-dependence 
of the statistics of 5Tr can be attributed to the vari- 
ations of Xt, and that the extension of Kolmogorov's 
RSH holds in the buoyancy-driven regime in turbulent 
Rayleigh-Benard convection. 

On the other hand, the Brandenburg shell model 
for homogeneous turbulent convection is, by construc- 
tion, free of boundaries and thus plumes. In an earlier 
work [30| , we have studied the anomalous scaling of an ac- 
tive scalar in the Brandenburg model and found that the 
scaling behavior is BO plus corrections due to the shcU- 
to-shell variations of the entropy transfer rate. These 
results thus verify the validity of the extension of Kraich- 
nan's RSH. We shall discuss these results in greater depth 
and also examine the validity of the extension of Kol- 
mogorov's RSH in Sec. E] 



IV. SHELL MODEL FOR HOMOGENEOUS 
AND ISOTROPIC TURBULENCE 

The basic idea of a shell model is to consider the veloc- 
ity variable u„, which can be associated with the veloc- 
ity difference Suj. with r = l/fc„, in discrete "shells" in 
Fourier space. Here fc„ = fcoA", with n = 0, 1, . . . , TV — 1, 
is the wavenumber of the nth shelljand A is customarily 
taken to be 2. In the Sabra model [5|, u„ is complex and 
satisfies the following equation of motion: 

at 



(19) 



where /(5„^o is the forcing acting only on the first shell, v 
is the kinematic viscosity and m* is the complex conjugate 
of Un- We use a = 1, 6 = —0.5, and c = —0.5. With this 



choice of the parameters, the model satisfies both energy 
and heUcity conservation in the inviscid Umit. 

Multiplying Eq. (fT9|) by u* , the complex conjugate of 
Un, and taking the average of the resulting equation and 
its complex conjugate, we obtain the energy budget for 
the nth shell: 



2 dt 
where 



Fu{kn) - F„(fc„+i) - J//cfju„P + "^(fuDSnfl , 

(20) 



Fu{kn) = fc„3«^iU>„+i + -u„<_i<_2) (21) 

Here SR and S represent the real and imaginary parts, 
respectively. The physical meaning of the different terms 
in the right hand side of Eq. (|20p is clear. The third term 
is the rate of energy dissipation in the nth shell due to 
viscosity, the last term is the power input due to external 
forcing, and F„ is the rate of energy transfer from the 
(n — l)th shell to the nth shell. In stationary state and 
in the inertial range at which the external forcing is not 
acting and energy dissipation is negligible, Eq. (|20|) gives 



« (F„(fc„)} - (F„(fc„+i)} 



(22) 



which implies that (-F'u(fc„)) is independent of fc„ in the 
inertial range, which is the statement of energy cascade. 
The ensemble averages are evaluated as time averages 
when the system is in the stationary state. In the shell 
model, the local energy transfer rate 11^ can thus be nat- 
urally identified as Fu{kn) with fc„ = 1/r: 



n^ 



Fuikn) 



(23) 



As for the local energy dissipation rate, it has to be 
defined accordingly in the shell model. From Eq. (pO)) . 
the total energy dissipation rate (at time t) in the shell 
model is 



N-l 



e{t) = ^ vkl\un\ 



(24) 



ri=0 



We define the analog of the local energy dissipation rate 
e„ in the shell model as 



N-l 



Cr ^ £r) 



E-fcn'l 



(25) 



again with fc„ = 1/r. With this definition, we have the 
nice result of e„(i) approaching e{t) in the limit of n ^ 
or equivalently in the limit of r -^ L. Summation of 
Eq. ([20]) from n > 1 to A^ - 1 gives 



1 d 
2di 



7V-1 



y^ |u^P = Fu{kn) - e„ , 



Thus in the stationary state, we have 
(Fuikn)) = (e„) 



(26) 



(27) 



but Fu{kn) 7^ e„ and, in general, the two quantities can 
have different statistical features. 

The velocity structure functions Qp{kn) are statistical 
moments of |u„| and Qp{kn) scales with kn with scaling 
exponents 7p: 



Qp(fc„) = {lunin ~ k--'^ 



(28) 



As reported in Ref. [5|] , the values of 7p are close to those 
obtained in experiments and deviate from the K41 val- 
ues of p/3. In this work, we are interested in examin- 
ing the validity of the RSHs for accounting this anoma- 
lous scaling behavior of the velocity structure functions. 
With r = l/fc„, the shell-model analog of Kolmogorov's 
RSH [Eq. ^] is written as: 



= ^'J'k-'^' 



(29) 



where is a random variable independent of n and statis- 
tically independent of e„. Equation ((29)l implies that the 
conditional velocity structure functions at fixed values of 
e„ would have simple K41 scaling: 



}p{kn) = (|u„|P |e„) 



-p/3 



(30) 



For simplicity of the notation, we suppress the depen- 
dence on e„ in the conditional velocity structure functions 
Qp. The notations for the other conditional structure 
functions will be simplified in the same fashion. Simi- 
larly, the shell-model analog of Kraichnan's RSH reads: 



Un = ^\Fu{kn)\'^''k-'^^ 



(31) 



where V' is a random variable independent of n and sta- 
tistically independent of Fu{kn). One implication of 
Eq. dni) is that (|M„pP) and {\Fu{kn)\P)k-P have the 
same scaling behavior. This implication has indeed been 
confirmed in the GOY [Slj and Sabra [S] models. From 
the definition of F„ [see Eq. (|211) ]. it is not surprising 
that \Fu{kn)\ would have the same scaling behavior as 
kn\un\'^, and thus the above implication could well be 
a direct consequence of the definition of F^ikn) in the 
shell model. However, we emphasize that Eq. ([3T|) is not 
merely a statement of the relation of the scaling expo- 
nents of (|Mn|'^^) and (|F„(fc„)|P). In particular, it has an- 
other important implication, namely, the conditional ve- 
locity structure functions at fixed values of Fu{kn) would 
have simple K41 scaling: 



Qlikn) = {\Un\P iFuikn)) ^ k~P/^ 



(32) 



Equation (|32p is nontrivial if u„ fluctuates even when 
Fuikn) is fixed at some given value, which would be the 
case when ^ is a random variable with a certain proba- 
bility distribution rather than a number of fixed value. 



V. SHELL MODEL FOR HOMOGENEOUS 
TURBULENT CONVECTION 

Homogeneous turbulent convection has been pro- 
posed [32] as a convective flow in a box, with periodic 



boundary conditions, driven by a constant temperature 
gradient along the vertical direction. In Boussinesq ap- 
proximation, the equations of motion read [33| : 



du ^ ^ 
— - + li • Vu 
at 

de 



ot 



Vp + i^V^a+agez (33) 

(34) 



with V • u = 0. Here, p is the pressure divided by the 
density, 6 ~ T — (To — /3z) is the deviation of temperature 
from a linear gradient —j3, Tq is the mean temperature of 
the fluid, and z is a unit vector in the vertical direction. 
A shell model for homogeneous turbulent convection 
driven by a temperature gradient was proposed by Bran- 
denburg 0. In Brandenburg's shell model, the velocity 
and temperature variables are real. We denote the veloc- 
ity variable as Vn, to distinguish it from u„ in the Sabra 
model, and the temperature variable as 9n- The equa- 
tions of motions are: 



dv„ 
~dt 

d0n 



vk^Vn 



— -r ^kj,. 



= agdn + Akn{vl^_^ - Aw„w„+i) 

+ Bkn{VnVn~l ~ )^v'^+^) (35) 

= /3u„ -I- A{kn{Vn-lOn-l " K+lVnOn+l) 

+ B{knVnOn~l ~ kn+iVn+l9n + l) (36) 



where A, B, A, and B are positive parameters. Earlier 
work showed that the scaling behavior depends only on 
the ratio B/A Q. It has been recently shown [30| that 
buoyancy drives the dynamics and affects the statistics of 
Vn and 9n when B/A is greater than some critical value of 
about 2. When buoyancy is driving the dynamics, energy 
is transferred from small to large scales on average [8| . 
As a result, a linear damping term —foVnSn,o [3J, |35| 
has to be added to Eq. ((35| for the system to achieve 
stationarity. 

In Bousinessq approximation, entropy is proportional 
to the volume integral of the temperature variance. En- 
tropy in the nth shell is, therefore, defined as 5„ = 0^/2. 
By studying the entropy budget obtained from Eq. ((36|) 
upon multiplication by 0„: 

-^ = Feikn) - Fg{k„+i) - Kklel + Pvnen (37) 

we get the rate of entropy transfer or entropy flux from 
{n — l)th to nth shell as: 



Fgikn) = kn{AVn-l + BVn)On- 



(38) 



In the stationary state and in the intermediate range in 
which (vnOn) and the entropy dissipation Kk'^{9'^) are 
neghgible, Eq. ([57]) gives: 



0^{Fe{kn))~{Fe{kn+i)) 



(39) 



Thus {Fg{kn)) is independent of A:„ in the intermediate 
range of scales, which is the statement of entropy cascade. 



We naturally identify the local entropy transfer rate Ilf 
as Fg{kn) with A;„ = l/r: 



n^ ^ Fg{kn) 



(40) 



In the Brandenburg shell model, the total entropy dissi- 
pation rate (at time t) is: 



x{t) = Y. nklel 

n=0 



(41) 



Thus as in Eq. ([^5]) . we define the analog of the local 
entropy dissipation rate in the shell model, %„, as 



Xr ^ x„ = ^ i^kl,el, 



(42) 



again with fc„ = l/r. 

The velocity and temperature structure functions, de- 
fined by 






{\dn?) 






(43) 
(44) 



have been studied recently (30| and found to have anoma- 
lous scaling behavior. We would like to examine whether 
this anomalous scaling can be understood using the 
RSHs extended to turbulent convection. The shell-model 
analogs of the RSH of Kolmogorov extended to turbulent 
convection, Eqs. (fTTI) and P^ . are: 



(45) 
(46) 



Vn = 


'^Mgf'xV'K^" 


On = 


'^eiagr^'xTK'" 


which imply that 




^p[kn) 


^ {\Vn\^ Xn)-k-'^^' 


Rp{kn) 


^ (i^«r xn)-k-p^' 



(47) 
(48) 



Similarly the shell model analogs of the RSH of Kraich- 
nan extended to turbulent convection, Eqs. pS]) and HI 
are 1301: 



^M9f^'\Feikn)\'/'k-y' 
^e{ag)-'^'\Fe{kn)\'/'k-'/' 



and they imply 



S;ikn) = {\Vnf |Fe(fc„))^fc-3p/5 
R;{kn) = {\en\P \Fg{kn)) ^ k-P/' 



(49) 

(50) 

(51) 
(52) 



Equations ([5T|l and (|52p have been confirmed in Rcf. f3(| 
thus supporting the validity of the extension of Kraich- 
nan's RSH to turbulent convection [Eqs. ((^ and ([50)1 ] 
in the Brandenburg model. The validity of the extension 
of Kolmogorov's RSH to turbulent convection [Eqs. flS)) 
and (|46p ] in the Brandenburg model will be examined 
and discussed in Sec. IVII 



VI. RESULTS AND DISCUSSIONS 

A. Validity of the RSHs of Kolmogorov and 
Kraichnan in Sabra model 

We numerically integrate Eq. (|19p using the fourth or- 

— 1/3 

der Runge-Kutta method. Starting with u„ — xkn , 
where x is a random variable, we evolve the equations in 
time for a short period of time. The results so obtained, 
with the phases randomized, are used as the initial data 
for the actual runs. Following Ref. 0, we take the exter- 
nal forcing / as a time correlated noise with exponential 
correlation and correlation time r such that it foUws the 
evolution equation: 

fit + dt) 
- f{t)E + aV-2(l - E^) logio(pi) exp(z2^p2) (53) 

where E = exp(— di/r), a is the standard deviation of /, 
and pi and p2 are two uniform random numbers between 
and 1. We use a = 0.01, r = 1, and f[t = 0) = 
5(1 + i) X 10-3. 




FIG. 1: The velocity structure functions Qpikn) as a function 
of the shell index n with p = 1 (diamonds), p = A (circles), 
p — 7 (squares), and p = 10 (triangles). The solid lines are 
the conditional velocity structure functions Qp(/c„) at given 
values of e,i with IQ-* < en < 3 x lO"* for the same values 
of p. 



We calculate the velocity structure functions Qp, and 
the conditional velocity structure functions Qp and Q* re- 
spectively at fixed values of e„ and i^„(^n)- The statistics 
are collected by averaging over a time of approximately 
1000 eddy turnover time of the largest scales. The scal- 
ing behavior of the model does not depend on the value 
of fcg and is also independent of the value of v when h' 
is small enough. We have used several different values of 



N. For the results reported below, we have used ko = 2, 
ly = 2x 10"® and N = 24, and i^ = 3.2 x 10"'^ and N = 21. 
Some of the results for the velocity structure functions Qp 
are shown in Fig. [1] We observe that they exhibit nice 
power-law dependence on fc„ for 3 < n < 16. The scaling 
exponents jp obtained are plotted in Fig. [51 These values, 
which are in good agreement with the values reported in 
the literature [5|, deviate fromp/3 demonstrating clearly 
that the velocity structure functions exhibit anomalous 
scaling. 
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FIG. 2: The scaling exponents yp (circles) for Qp{kn) as a 
function of p. The solid line is the K41 result of p/3. The 
error increases with p and we show the largest error for p = 10. 

Our aim is to examine the validity of Kolmogorov and 
Kraichnan's RSHs [Eqs. ((29l) and ((3T|) ] in accounting for 
the origin of the anomalous scaling behavior by examin- 
ing the validity of Eqs. pO)) and (|32p. To do so, we study 
the scaling exponents of Qp{k„) and Qp{k„), defined by: 



Qp\kn) 



n 



(54) 
(55) 



and check directly whether 7p and 7* agree with p/3. 

To calculate the conditional velocity structure func- 
tions Qp at fixed values of £«, we average only those |un|^ 
when the value of e„ is within a given narrow range of 
values. The results for Qp are also shown in Fig. [T] It 
can be seen that Qp are different from Qp. Moreover, the 
scaling range of Qp is shorter than that of Qp. Thus |m„| 
is correlated with e„. Yet the the scaling exponents 7^ 
and 7p are close to one another. In particular, as shown 
in Fig.m 7^ continue to deviate from p/3, demonstrating 
that Eq. j^ and thus Eq. ^ is invalid. 



Similarly, to calculate the conditional velocity struc- 
ture functions Q* at fixed values of Fu(kn), we average 




FIG. 3: The scaling exponents 7p (triangles) an 7p (circles) of 
Qp{kn) and Qp{kn) respectively as a function of p. The solid 
line is the K41 result of p/3 and the largest errors for p = 10 
are shown. 
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FIG. 5: Compensated plots of Qp{k,i)kn versus n. Same 
symbols as in Fig. U) 



only those \un\^ when Fu{kn) assumes values in the same 
given narrow range. Some of the results are shown in 
Fig. m It can be seen that Q* are again different from Qp 
thus showing that m„ and Fu{kn) are correlated. More- 
over, the scaling range of Qp is comparable to that of Qp. 
In Fig. [H we see that the scaling exponents 7* are consis- 
tent with p/3. We also present the compensated plots of 

Qp{kn)kn in Fig. [5] It can be seen that the data points 

are consistent with being independent of n in the scaling 

range, thus demonstrating directly the scaling behavior 

73/3 

of Q* ^ kn ■ These results thus confirm Eq. 



density functions of u„ at fixed values oi Fu{kn) and con- 
firm that the distributions indeed have a finite extent [see 
Fig. [B]. Moreover, the simple scaling behavior of Q* im- 
plies that the conditional probability density functions 

P{Yn\Fu{kn)), with y„ EE |u„ |/y^(M2 |F„(fc„)) , are scale 
invariant or n- independent. This is indeed the case as 
seen in Fig. [SI Our results thus show that Kraichnan's 
RSH [Eq. dnil)] is valid in the Sabra shell model. 





FIG. 4: The conditional velocity structure functions Qp(fcn) 
at fixed values of F„(fc„) with 10"* < -F„(fc„) < 3 x IQ-". 
Same symbols as in Fig. [T] 



FIG. 6: The conditional probability density functions 
PiYn\Fu{kn)) for n = 13 (circles), n = 15 (squares), 
n = 17 (triangles). 



As discussed in Sec. IIVI it is important to check that 
Un indeed fluctuates even when Fu(kn) is fixed at some 
given values. Thus we study the conditional probability 



B. Validity of extensions of Kolmogorov and 

Kraichnan's RSHs to turbulent convection in 

Brandenburg model 

We numerically integrate Eqs. ([55]) and ([M]) using 
fourth-order Runge-Kutta method with an initial con- 
dition of Vn = 9n = except for a small perturbation 
of 6n in an intermediate value of n. We use ko — 1, 
A = 0.01,5 = 1, /? = 1, i = B = 1, ag = 1, 
1/ = 5 X IQ-i^ AC = 5 X 10-15, /o = 0.5, and N = 32. We 
calculate the velocity and temperature structure func- 
tions Sp and Rp when the system is in stationary state. 
The results for 5*^ and Rp are shown in Figs. [7| and [5] re- 
spectively. It can be seen that they all have good scaling 
behavior. 
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FIG. 7: Sp for p = 1 (circles), p = 3 (squares), p = 6 (trian- 
gles) 
for the same values of p. 



The solid and dashed lines are respectively Sp and Sp 



FIG. 8: Rp for p = 2 (circles), p = 5 (squares), p = 8 (trian- 
gles) . The solid lines are Rp for the same values of p. 



0.2 



lo 0.( 



djOSOQgf 



-0.2 



JJLP 



id -0.4 

m 



-0.6 



-0.8. 



Jn "'"''^^OOr^ 



D 



D 



D 



n 



D 



D 



D 



D 



[] 



4 6 8 10 

P 



FIG. 9: The deviations of C,p and ^p from the BO values: 
Cp — 3p/5 (circles) and ^p — p/5 (squares). 



In an earlier study [30| , we have studied the scaling be- 
havior of an active scalar using the Brandenburg model 
and found that C,p and ^p are given by the BO values plus 
corrections (see Figl^]). This shows that 5^ and Rp have 
anomalous scaling behavior. In the same study, we have 
checked directly that the scaling exponents of S* and R* 
denoted by Cp and ^*, have the BO values of 3p/5 and 
p/5, respectively, thus confirming the validity of Eqs. (j5ip 
and (152]) . To calculate S* and R* we average \vnY' a-nd 
|6I„|P only when Fe{kn) = 0.2 ± lO^^. We show S*p also 
in Fig. [71 It can be seen that S* are close to Sp but with 
a slightly steeper decrease with fc„. Results for R* are 
shown in Fig. [in] in which the steeper decrease with kn 
can be seen more clearly. The validity of Eqs. (jSTj) and 
([5^ indicates the validity of the extensions of Kraich- 
nan's RSH to turbulent convection in the Brandenburg 
model. 



With these results, we expect that the extensions of 
Kolmogorov's RSH to turbulent convection would be in- 
valid in the Brandenburg model. In particular, we expect 
Eqs. (|i7|) and (|i5)) to be invalid. To show this, we cal- 
culate Sp and Rp by averaging \vnY' and \9n\^ only when 
Xn = 0.2 ± 10"^. Results for Sp and Rp are shown also 
in Figs. [7] and [SI respectively. We find that both Sp and 
Rp are very close to Sp and Rp with Rp having a shorter 
scaling range than Rp. The scaling exponents of Sp, de- 
noted by C,p, are close to Cp which are close to 3p/5 within 
errors. So we focus on the scaling exponents ^p of Rp. In 
Fig. [Til we compare S,p with ^* and the BO values of p/5 
and see clearly that, as expected, Eq. (|48l) does not hold. 
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FIG. 10: The conditional temperature structure functions Rp 
at fixed values of Feikn) for p — 2 (circles), p = 5 (squares), 
p = 8 (triangles) . 
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FIG. 11: Comparison of ^p (triangles) and Q (circles) with 
the BO values of p/5 (solid line). The largest errors for p = 10 
are shown. 



VII. CONCLUSIONS 

One longstanding goal in turbulence research is to un- 
derstand, from first principles, the origin of the anoma- 
lous scaling of the fluctuating physical quantities. An im- 
portant idea was proposed by Kolmogorov in his refined 
theory [3| , which attributes the origin of the anomalous 
scaling behavior of the velocity fluctuations in homoge- 
neous and isotropic turbulence to variations of the lo- 
cal energy dissipation rate. It was pointed out later by 
Kraichnan [3| that it is more appropriate to replace the 
local energy dissipation rate by the local energy transfer 
rate. We refer these two ideas to as Kolmogorov's and 
Kraichnan's RSHs. Both of them can be extended to 
account for the anomalous scaling behavior of the veloc- 
ity and temperature fluctuations in turbulent convection 
in which temperature acts as an active scalar. Specifi- 



cally, in the extension of Kolmogorov's RSH to turbulent 
convection, the local entropy (or temperature variance) 
dissipation rate plays the role of the local energy dissipa- 
tion rate. Similarly, in the extension of Kraichnan's RSH 
to turbulent convection, the local entropy transfer rate 
plays the part of the local energy transfer rate. 

In this paper, we have examined the validity of Kol- 
mogorov's and Kraichnan's RSHs and their extensions 
to turbulent convection respectively in the Sabra shell 
model of homogeneous and isotropic turbulence W] and 
in the Brandenburg shell model of homogeneous turbu- 
lent convection @]. The validity of Kolmogorov's RSH 
has been examined in previous studies. These studies 
mainly focussed on the statistical correlation of the ve- 
locity difference and the local energy dissipation rate and 
the relation between the scaling exponents of the velocity 
structure functions and the moments of the local energy 
dissipation. The validity of Kraichnan's RSH is much less 
studied and again the focus is on the relation between the 
scaling exponents of the velocity structure functions and 
the moments of the local energy transfer rate. An impor- 
tant consequence of Kolmogorov's or Kraichnan's RSH is 
that the conditional velocity structure functions at given 
values of the local energy dissipation rate or the local 
energy transfer rate would have simple K41 scaling be- 
havior. Similarly, an important consequence of the exten- 
sions of Kolmogorov's or Kraichnan's RSH to turbulent 
convection is that the conditional velocity and temper- 
ature structure functions at given local entropy dissipa- 
tion rate or the local entropy transfer rate would have 
simple BO scaling behavior. Our approach is to check 
directly these consequences. In the shell models, the lo- 
cal energy or entropy transfer rate is easily identified with 
the shell-to-shell energy or entropy transfer rate. On the 
other hand, the local energy or entropy dissipation rate 
has to be defined accordingly. As shown in Fig. [3l the 
scaling exponents 7* of the conditional velocity struc- 
ture functions at given shell-to-shell energy transfer rate 
indeed have the K41 values of p/3 while the scaling expo- 
nents 7p of the conditional velocity structure functions at 
given shell-model analog of the local energy dissipation 
rate continue to deviate from the K41 values. This re- 
sult shows that Kraichnan's RSH but not Kolmogorov's 
RSH holds in the Sabra model. Similarly, as shown in 
Fig. [Ill we have found that the scaling exponents ^* of 
the conditional temperature structure functions at given 
shell-to-shell entropy transfer rate have the BO values of 
p/5 while the scaling exponents ^p of the conditional tem- 
perature structure functions at given shell-model analog 
of the local entropy dissipation rate continue to deviate 
from the BO values. Our result thus shows that the ex- 
tension of Kraichnan's RSH but not Kolmogorov's RSH 
to turbulent covection holds in the Brandenburg model. 
In summary, our work shows that Kraichnan's RSH and 
its extension to turbulent convection hold in shell models 
of turbulence. It would be interesting to perform similar 
study in direct numerical simulations and in experimen- 
tal investigations. 
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